///'Vb 

NASA Technical Memorandum 106891 _ syyye % 

WP 


Electromagnetic Radiation in the Plasma 
Environment Around the Shuttle 


Boris V. Vayner and Dale C. Ferguson 

Lewis Research Center _ 

Cleveland, Ohio „ 

( NAS A-TM— 1 0689 1 ) ELECTROMAGNETIC N95-23116 ^ 

RA0IATI0N IN THE PLASMA ENVIRONMENT 

AROUND THE SHUTTLE (NASA. Lewis :1 

Research Center) 21 p Unclas 


&3#l/46 0044403 


Prepared for the . 

Chapman Conferenc e on Measurement Techn iques for Space Plasma 
sponsored by the American Geophysical Union 
Sante Fe, New Mexico, April 3-7, 1995 



National Aeronautics and 
Space Administration 



ELECTROMAGNETIC RADIATION IN THE PLASMA 
ENVIRONMENT AROUND THE SHUTTLE 


Boris V. Vayner*, and Dale C. Ferguson* 

National Aeronautics and Space Administration 
Lewis Research Center, Cleveland, Ohio 44135 

Abstract - 

As a part of the SAMPIE (The Solar Array Module Plasma Interaction Experiment) 
program, the Langmuir probe (LP) was employed to measure plasma characteristics 
during the flight STS-62. The whole set of data could be divided into two parts : I) low 
frequency sweeps to determine voltage-current characteristics and to find electron 
temperature and number density; ii) high frequency turbulence (HFT dwells) data caused 
by electromagnetic noise around the shuttle. The broadband noise was observed at 
frequencies 250-20,000 Hz. Measurements were performed in ram conditions; thus, it 
seems reasonable to believe that the influence of spacecraft operations on plasma 
parameters was minimized. The average spectrum of fluctuations is in agreement with 
theoretical predictions. According to purposes of SAMPIE, the samples of solar cells 
were placed in the cargo bay of the shuttle, and high negative bias voltages were applied 
to them to initiate arcing between these cells and surrounding plasma. The arcing onset 
was registered by special counters, and data were obtained that included the amplitudes of 
current, duration of each arc, and the number of arcs per one experiment. The LP data 
were analyzed for two different situations: with arcing and without arcing. Electrostatic 
noise spectra for both situations and theoretical explanation of the observed features are 
presented in this report. 

1. Introduction 

During the space shuttle flight STS-62 a broadband electrostatic noise was observed at 
frequencies 250-20,000 Hz. The measurements were performed in ram conditions by the 
Langmuir probe (LP) flown as a part of the SAMPIE package ( Ferguson et aL, 1994). 
The noise has almost a flat spectrum with a sharp decline near the lower hybrid resonance 
frequency / jj,, and the intensities are ranging from 0.1 to 5 mV/m ( Vayner & Ferguson, 
1995). It should be noted that such a noise was observed in ram conditions more than ten 
years ago (Siskind et aL, 1984), but there is not a satisfactory explanation for mechanism 
of generation yet. An investigation of the plasma wave turbulence within the wide range of 
frequencies (from a few Hz to 200 kHz) was done by using the Plasma Diagnostic 
Package (PDP) during the SL-2 and STS-3 flights (Murphy et aL, 1986; Gurnett et aL, 
1988) . It was established that the noise is electrostatic, the highest intensities occured in 
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the region downstream of the spacecraft, and the intensities increased considerably after 
water dumps {Pickett et al, 1989). Two different processes were considered in attempting 
to understand the generation (instability) of the plasma turbulence caused by water release: 
the drift instability and the Ott-Farley instability caused by the ring distribution of the 
water ions. Model computations of the ion distribution function confirm the idea that this 
function is non-maxwellian {Paterson & Frank, 1989). Moreover, so called “ring 
distribution” was measured directly, and one may believe that the electrostatic broadband 
noise observed in the wake of the shuttle can be explained theoretically {Kurth & Frank, 
199(f). However, an additional analysis of electrostatic noise generation mechanism is 
needed because there are essential differences in plasma parameters between SAMPIE 
and PDP experiments. As it was mentioned above that all the measurements of 
electrostatic noise (HFT dwells) were performed in ram conditions. The LP uses a 5 cm 
diameter spherical sensor mounted on a fixed boom approximately 100 cm from the 
surface of the shuttle {Hillard &. Ferguson, 1993; Morton et al, 1995), and fluctuations of 
electric current were measured by using the logarithmic amplifier with V-A 
characteristic shown in Fig. 1 {Bozich, 1994). Each dwell lasted for 4 ms, and 
experimental data were represented in digital form: 80 measurements of current with 50 
ps intervals. 

The spectra of fluctuations were computed by using the FFT procedure; thus, we are 
able to study the plasma turbulence within the narrow frequency interval Af=250-10,000 
Hz. According to the purposes of the SAMPIE, the samples of solar cells were placed in 
the cargo bay of the shuttle, and high negative bias voltages were applied to them to 
initiate arcing between these cells and surrounding plasma. It was shown that there is 
influence of arcing on the parameters of the electrostatic noise ( Vayner & Ferguson, 
1995), and we discuss this topic in the last chapter of the current paper. 


2. Observations 


Three examples of HFT dwells with their spectra are shown in Fig. 2-4. The 
background plasma parameters (electron number density and temperature) were obtained 
by using the V-A characteristic of the LP {Morton et al., 1995). As can be seen, there are 
no essential differencies among these graphs besides the amplitudes of fluctuations: the 
amplitude of the current reaches its magnitude i^ - 0.1 pA for dwell shown in Fig. 2, and 
the amplitudes are substantially higher for dwells shown in Fig.3 and Fig.4 {i^ » 0.4 pA). 
But if we take into account difference in electron number densities for these three dwells, 

6/i 

we may conclude that the level of fluctuations — — is five times higher for the last dwell 

n c 

(Fig.4). 

All spectra demonstrate sharp decline for frequencies fe 6 kHz. This fact could be 
considered as an argument in favor of the Ion Acoustic Waves hypothesis because such 
kind of waves can be excited within the frequency range F-<£< if the electron 
temperature T e is much greater than the ion temperature T- t . {F- t is the ion 
gyrofrequency). HFT dwell that was recorded one minute before the dwell shown in 
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Fig.3 demonstrates almost the same signal shape and spectrum but the amplitude of 
fluctuations is about five times less (Fig.5). 

One minute (60 seconds) is not a natural time scale for IAW because the calculations of 
the attenuation rate and instability increment for such waves show very short time intervals 
(x=3-100 ms) (Fig.6). The duration of one dwell (td= 4 ms) is not long enough to observe 
any regular trend in amplitudes although one could believe that the decrease of the 
am plitude for time interval t>2 ms is caused by the real attenuation (Fig. 3 and Fig. 4). 
Some measurements were done in plasma with relatively low electron number density 
(Fig.7). The pressure of neutral gas is about 20% less than for the dwells shown above 
but surprisingly the level of fluctuations is almost one order in magnitude higher. There are 
no essential differences in spectra between this dwell and other dwells . 

For very low electron number density and low presssure of the neutral gas the 
instrument noise was registered only (Fig.8). The origin and the characteristics of the 
observed electrostatic noise are needed in theoretical explanation. We interpret these 
fluctuations in terms of the Ion Acoustic Waves travelling in low ionized plasma layer 
surrounding the shuttle. 

3. Interpretation 


The shuttle is surrounded by the gas cloud caused by many different processes 
accompanying the spacecraft operation: outgasssing of surfaces, leakages in valves, 
thrusters firing, etc. (Paterson & Frank, 1989). Molecules of the gas are moving at the 
thermal speed within the rest frame of the shuttle, and ionospheric ions and neutral atoms 
can be considered an inflow with an average speed Vs=7.8 km/s. If we suggest that the 
cloud is formed from water vapor (in the main part) with initial temperature 
T s 300 K , we may determine the number density of water molecules from 


measurements of the gas pressure near the shuttle: 



( 1 ) 


This value of the number density is few times more than the number density of the 
ambient atmospheric gas at flight altitudes (220-310 km). According to direct 
computations, the water vapor densities at distances about 50 m from the shuttle can be 
as large as n(/J 2 o) - 2 • 10 9 cm' 3 , and this value reaches 4 -10 10 cm' 3 in the 
nearest vicinity of the spacecraft surface (Paterson & Frank, 1989). Now, we can 
estimate the mean free path for the water molecule. The thermal speed is equal to 
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If we adopt the cross section for molecular collisions a - 2- 10 -15 cm 2 ( Kaplan & 

PikeVner, 1979) and the number density of the neutral gas n g ■* 10 10 cm' 3 , we will get a 
rough estimation for the mean free path: 

/ - 5- 10* cm - 05 km (3) 

O'l t g 

that is much greater than the linear dimension of the spacecraft (L^-10 m). The time 
interval between two collisions is equal to 



These estimations allow us to write the radial dependence of the water molecules number 
density in the simple form ( Paterson & Frank, 1989): 

(5) 

The measurements of the electron temperature in the vicinity of the shuttle show the 
magnitudes inside the narrow interval T e =(0. 1-0.3) eV ( Morton et aL, 199S). For the 
IAW, the phase velosity is equal to 



According to the dispersion relation for IAW, the wave length can be calculated as 



where/ is a frequency of fluctuations measured in the rest frame of the plasma. 

Thus, within the framework of IAW hypothesis the wave length of fluctuations is less than 
3 m for whole range of the observed frequencies. It means that we may consider the gas 
layer surrounding the shuttle as slightly non-homogenious because of the following 
inequalities: 

Su>*K» <<L<<l ( 8 ) 

In the relations (8) Su> is equal to the distance of LP from the spacecraft surface. 

Neutral water molecules that are ejected from the shuttle surfaces will be ionized by an 
electrical charge exchange with ions 0 + having the flow velocity Vo=- 7.8 km/s and 
temperature T+slOOO K: 

O* +H 2 0-+H 2 0* +0 (9) 
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Three more reactions should be added to the reaction (9) to write the chain of reactions 


correctly: 

H 2 0 + + H 2 0 H 3 0 * + OH (a) 

H 2 0 + + e’ — H 2 0' -+H + OH ( b ) (10) 

H 3 0* +e~ -* H + H 2 0 and OH + 2 H (c) 

The kinetic equations for the reactions (9) and (10) were solved by using reaction rates 
tabulated in the work of Paterson and Frank (1989). For particular case 
n(o+) - 2- 10 5 cm' 3 , 

n(H 2 o) - 6 • 10 10 cm 3 ,and n € - n(o * ) + n(H 2 0 * ) + n(H 3 0+ ) 
we obtain the ratio 


X m 


AW) 

Aw) 


3.9 • 10" 6 


( 11 ) 


As a result, the plasma comoving the shuttle can be characterized by the ion number 
density as large as (2 - 3) * 10 s cm' 3 . In the rest frame of the shuttle these ions are 
undergone of the influence of magnetic field B and induced electric field E that is 
perpendicular to both vectors V, and B, and the electric field strength is equal to 


E -V, -21 -sin© - 0.23- 



•sin© 


V 
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( 12 ) 


where 0 is the angle between the Earth magnetic field and the velocity of the spacecraft 
(Fig.9). 

It is well known that electrically charged particles will drift with speed V D - V t ■ sin© in 
the direction that is perpendicular to both electric and magnetic fields. The angle between 
the shuttle velocity and the magnetic field had been varying from ©=-65° to ©=90° during 
the experiments that are considered in the current paper. It means that the projection of 
drift velocity directed to the surface of the shuttle V± is always greater than the 
projection along the surface V\ \ : 


^ L -|tan©|>l (13) 

*// 

We can conclude that there are two streams of ions with relative velocity 


Vr=-V t - V D . (14) 

Thus, it is possible to use the dispersion equation for IAW that may be written in the 
following form (Akhiezer et al, 1973 ): 
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where A 2 - — M — is the Debye length, co^ are the ion plasma frequencies for 

An 'e n c 

oxygen and water ions respectively, and a is the angle between the wave vector k and 
the relative velocity V, . 

The fourth order equation (15) has complex roots (instability) if the following inequality 
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can be fulfilled: 


Vj cos 2 a <Wj 


>» W +<] 

CD 1 +0) 2 


To obtain the inequality (16) we suggest T tl - T e2 , and we take into account that masses 
of oxygen and water ions are almost equal to each other. 

For the simplest example (o x « o > 2 , the condition of instability (16) can be written in the 

form 

W 

|cosa|<2~, (17) 

V 


where V=l^cos 0 . 

It is seen from the expression (17) that the cone of the instability is wide for the angle 0 
near the right angle. But the dispersion relation does not have complex roots when angle 
©=ji/ 2. It is obvious because in this case the drift velocity is equal to the space shuttle 
speed exactly, and there is no any relative motion of two kinds of ions; thus, the two- 
stream instability does not work when the spacecraft velocity is directed perpendicularly to 
the Earth magnetic field. 

Now we can consider some particular experiments. For example, angle between vectors 
V r and B is equal to 0= - 66.4° (Fig.2). It is easy to calculate the IAW phase speed 
Wr=8-10 4 cm/s and |cosa | < 05 . The complex frequency (the growing mode only) can be 
determined from the Eq.(15): 


- 05+ i- 034 (cosa - 035) 
kW t 

- 77 - - 0.9 + i • 034 (cosa - 0.45) 

kW t 


(18) 


U>1 

For different ion number densities — - 2 we obtain almost the same result: 

to 2 

■^ 7 - - 0.6 + 1 • 033 (cosa - 035) (19) 

kW s 

During the experiment E_62-2/03 the electron temperature was substantially higher 
(Fig.3). IAW speed was equal to W' x =9.5-10 4 cm/s, and the relative velocity of the ion 
streams was equal to V=1.36-10 5 cm/s. It is seen that for this particular experiment 
“beam” velocity is almost equal (a little higher) to the IAW phase velocity. Such 
condition is optimal for the energy transfer from flux kinetic energy to fluctuations, and 
the rate of instability is high for all the magnitudes of angle between wave vector and 
stream velocity (see Eq. (17)). The result depends on the relation between ion number 
densities rather slightly: 
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When the angle between the spacecraft velocity and the magnetic field is almost equal to 
90° the relative speed is less than the IAW phase speed, and the conditions for the two- 
stream instability are broken in this case. The dispersion relation (15) has two real roots 


for any values o>i and C 02 : 
co 


kW. 


±1 


( 21 ) 


Due to the IAW attenuation, we should not observe the electrostatic noise during the 
time interval t— 5 min which is needed to change the angle 0 from 85 to 95 degrees 
(Fig. 10). However, fluctuations were observed even during the dwell recorded for V,1B 
(Fig.4). Such situation was discussed by Pickett et al.(1989) , and they suggested that a 
modified two-stream instability (MTSI) could explain the observational data obtained in 
their experiment for the wake of the shuttle. MTSI generates fluctuations with 
frequencies near the lower hybrid resonance frequency (few kHz), and some nonlinear 
processes should be considered to explain the observed spectrum that has significant 
amplitudes at frequencies as low as a few Hz. 

If we confront our LP data obtained for Br = 0 (Flg.4) with data obtained for Br&\B\ 
(Fig.3 and Fig.5), we may su gg est that the nature of the turbulence is the same in all three 
cases because the shape of signals and their spectra do not have any considerable 
differences. One of the possibilities to explain the IAW instability is the gradient instability 
caused by electron and ion number densities gradient with a scale 


1 n ; 


dn i \ 1 

dR, 


>5m (see Eq. (5)). 


following form (Timofeev, 1964): 


The dispersion relation can be written in the 
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where (p is the angle between the wave vector k and the Earth magnetic field B. 


Taking into account that the Debye length is much smaller than the wave length, we can 
represent the increment of instability in the following form: 


Y_ 

to 
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( 23 ) 


where <0 obeys the Eq. (21). 
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It should be noted that the Eq.(23) is valid when the angle between the wave vector k and 

k k 

the magnetic field B is not close to zero, but — < 1 ,and ~ <1 ( ky is the projection of 


the wave vector on the direction that is perpendicular to both vectors B and grad n,.). 
The ratio (23) can be easily estimated for the considered frequency range: yAo= 0.005-0.1. 
The attenuation rate of IAW depends on the ratio of the electron temperature to the ion 
temperature ( Krall& Trivelpiece, 1973) (Fig. 11). The electron temperature is measured 
for each dwell, and we can use these data for further calculations. 

We can estimate the ion temperature working within the framework of our hypothesis 
that the main ion component of the plasma is water ions generated by charge exchange 
process. Thus, we can believe that T-, <300 K. Moreover, the IAW attenuation rate 
depends on the angle between the wave vector and the magnetic field ( Alpert, 1990) 
(Fig.6). If we compare attenuation rates mentioned above with increment (23), we see 
that there is enough room for instability of IAW even in ram conditions with V S 1B if 



4. The Influence of Arcing 

As it was mentioned above HFT data were obtained during the experiments when arcs 
occured. First of all, it is seen that there is a significant difference in amplitudes between 
signals with arcing and without them (Fig. 12). For this particular dwell, the amplitude of 
current fluctuations is as large as 2 ^A. The number of arcs, the maximum current for 
each arc, and the duration of arc were measured by the special counter. One example of 
data obtained during the experiment E_60-l/01 is shown in Fig. 13. It was established 
that there is a strong correlation between the number of arcs observed in an experiment 
and the amplitude of fluctuations (Fig. 14). All these facts show that there is an influence 
of arcs on HFT data. We believe that the large amplitude of LP current is caused by 
changes in the spacecraft potential due to arcing ( Vayner & Ferguson, 1995). However, 
this problem needs in more profound investigation. The measurements of the time 
interval between arc and HFT dwell would be very useful in this respect. 

5. Conclusion 

During the Solar Array Module Plasma Interactions Experiment (SAMPIE) fluctuations 
in the electron number density were observed in ram conditions. The amplitudes of these 
fluctuations were as large as 10' 3 - KT 4 from the background density at frequencies 
f= 250-6500 Hz. The measurements of the Earth magnetic field allow us to show that 
this frequency range is inside the interval F i < f < f w ; thus, we might suggest that the 
Ion Acoustic Waves were observed. The additional argument in favor of the hypothesis of 
IAW is the inequality between ion and electron temperatures. From having measurements 
there was determined that the electron temperature was varying within the interval 
O.l^r^O.3 eV {Morton et al,1995). According to Patterson and Frank (1989) the ion 
temperature should be substantially lower : 7i < 0.03 eV. This value of ion temperature 
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is in accordance with the measurement of the pressure of neutral gas around the shuttle 
(Fig.15) and our estimations of the number density of water molecules (see also Eq.(l)). 
The facts nam ed above allow us to adopt the basic model of the gas layer around the 
shuttle that was elaborated by Paterson and Frank ( 1989 ). Two kinds of instabilities were 
su gg ested for explanation of the origin of IAW: two-stream instability that works when 
the angle between the spacecraft velocity and the magnetic field is not too close to a right 
angle, and the gradient instability that works for situation when the magnetic field and 
velocity of the shuttle are perpendicular to each other. It should be noted that the last 
instability is weaker than the first one, and we could observe the transition between the 
two regimes if the dwell would last for several minutes at least More measurements need 
to be done, particularly for frequency range 20-100,000 Hz, for understanding of the 
nature of the observed turbulence. 
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Fig.l Voltage-current characteristic of the Logarithmic Amplifier. 
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Fig. 2a HFT dwell recorded during the experiment E_44-2/01. MET 7/15:43, 
electron number density 5.4- 10 s cm' 3 , electron temperature 0.12 eV, 
and neutral gas pressure 2.27 microTorr. 
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Fig.3a Experiment E_62-2/03. MET 7/23:40, n e =4.7*10 5 cm* 3 , T e =0.17 eV, and 
p=2.1 pTorr. 
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Fig. 4a Experiment E_48-2/01. MET 7/22:54, n e =6 10 s cm 3 , T e =0.1 eV, and 
p=1.86 pTorr. 
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Fig. 5a Experiment 62_2/02. MET 7/23:39, n e =3.8 - 10 s cm' 3 , T e =0.17 eV, and 
p=2.2 (xTorr. All parameters are almost the same as in Fig. 3, but the 
amplitude of fluctuations is about five times less. 
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Fig. 5b. This spectrum demonstrates very sharp decline for frequencies hfiR. 



Fig. 6 IAW attenuation time vs. frequency for three magnitudes of angle between 
the Earth magnetic field and the wave vector (<J>). (T,/T|=15) 
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p=1.85 |iTorr. 



Fig. 7b The spectrum of the signal shown in Fig. 7a. B=033 Gs, and 0=74°. 
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LP VOLTAGE (V) 



Fig. 8 Experiment E_52-12. MET 10/06:27, n.=310 4 cm 3 , T,=0.1 eV, and 
p= 0.665 pTorr. The instrument noise is registered for low electron 
number density and low pressure of the neutral gas. 



Fig. 9 The spacecraft velocity is directed along x axis. The vector of magnetic field 
is placed on xOz plane, and the induced electrical field is perpendicular to 
both vectors B and V,. The vector of drift velocity V D is perpendicular to 
both vectors B and E. 
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Fig.12 LP output voltage fluctuations measured during the dwell in ram conditions 
when arcing occurred (MET 7/17:21). This plot demonstrates the large 
amplitude caused by arcing. 



Fig. 13 Data from the Experiment #60-1/01. The bias voltage V=-400 V , and 181 
arcs were registered. 
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Fig. 14. LP peak voltage vs. the number of arcs per one experiment Twelve dwells 
are used to draw this chart No correlations were found between the 
peak voltage and the distance of the LP from biased solar cell samples. 
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Fig. 15 The neutral gas pressure vs. time started at MET 7/22:23 (SAMPIE) 
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reasonable to believe that the influence of spacecraft operations on plasma parameters was minimized. The average 
spectrum of fluctuations is in agreement with theoretical predictions. According to purposes of SAMPIE, the samples of 
solar cells were placed in the cargo bay of the shuttle, and high negative bias voltages were applied to them to initiate 
arcing between these cells and surrounding plasma. The arcing onset was registered by special counters, and data were 
obtained that included the amplitudes of current, duration of each arc, and the number of arcs per one experiment The LP 
data were analyzed for two different situations: with arcing and without arcing. Electrostatic noise spectra for both 
situations and theoretical explanation of the observed features are presented in this report. 
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